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P^ . ABSTRACT 

Some years ago Dray and 't Hooft found the necessary and sufficient conditions to 
Oh! 

^ . introduce a gravitational shock wave in a particular class of vacuum solutions to Einstein's 



equations. We extend this work to cover cases where non-vanishing matter fields and 
cosmological constant are present. The sources of gravitational waves are massless particles 
moving along a null surface such as a horizon in the case of black holes. After we discuss the 
general case we give many explicit examples. Among them are the (i- dimensional charged 
black hole (that includes the 4-dimensional Reissner-Nordstrom and the (i- dimensional 
Schwarzschild solution as subcases), the 4-dimensional De-Sitter and Anti-De-Sitter spaces 
(and the Schwarzschild-De-Sitter black hole), the 3-dimensional Anti-De-Sitter black hole, 
as well as backgrounds with a covariantly constant null Killing vector. We also address 
the analogous problem for string inspired gravitational solutions and give a few examples. 
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1. Introduction 

In an S'-matrix approach to black hole physics one has to take into account the 
interactions between Hawking particles emitted from the black hole. However in order 
to do that one would like also to know how the presence of particles around a black hole 
affects the black hole itself (the black hole reacts back). It has been convincingly argued 
that the gravitational interaction between particles at Planckian energies (see 
dominates any other type of interaction and therefore one needs to know what gravitational 
effects the particles have on the original black hole geometry. 

For the case of a massless particle moving along the horizon of a Schwarzschild black 
hole the result of this elementary classical black hole back reaction has been found by 
Dray and 't Hooft . What they found generalized the gravitational shock wave due to a 
massless particle moving in fiat Minkowski space [0 . In fact contains the necessary and 
sufficient conditions for being able to introduce a gravitational shock wave via a coordinate 
shift in more general vacuum solutions to Einstein's equations.0 

The purpose of this paper is to extend the results of to the case when matter 
fields and even a non-vanishing cosmological constant are present since, many interesting 
gravitational solutions belong in this category. The organization of this paper and some of 
the important results are: In section 2 we address the problem of how a quite general class 
of solutions to Einstein's equations in the presence of matter and cosmological constant 
(in d dimensions) back reacts to a massless particle moving along a null surface. We also 
examine the geodesies of particles moving in these geometries and show that discontinuity 
and refraction effects take place as they try to cross the null surface. Our treatment, as 



^ The solution of can be formulated as two Schwarzschild black holes of equal masses glued 
together at the horizon. For a spherical shell of massless matter moving along u = uq ^ Q the 
solution [^] represents two Schwarzschild black holes of unequal masses glued together at u = uq. 
Various gravitational shock waves obtained by infinitely boosting Q @ known solutions have been 



found (see for instance g @ li Q ) . 



the one in |^], is non-perturbative but exact. In section 3 we consider the case of the d- 
dimensional Reissner-Nordstrom charged black hole fl^. After we give the general result 



we concentrate to the 4-dimensional case which physically is the most interesting one. We 
find that as we approach the extremal case of equal mass and charge the effect the massless 
particle has become gigantic in magnitude. In section 4 we consider firstly the De-Sitter 
and Anti-De-Sitter constant curvature spaces in four dimensions. In the former case we 
find that the discontinuity in geodesies we have already mentioned depends in a rather 
unexpected way on the angular distance from the position of the massless particle. We 
also consider the case of the Schwarzschild-De-Sitter black hole which interpolates between 
the De-Sitter space and the Schwarzschild black hole. We close this section with the 3- 
dimensional Anti-De-Sitter space which, after a discrete identification, can be interpreted 



as a 3-dimensional black hole [|T3[. In section 5 we consider the analogous problem for 
string inspired gravitational solutions. In particular for a 4-dimensional electrically and 
magnetically charged dilatonic black hole [|14|[^ and the background corresponding to 



the conformal field theory SL{2, IR) /IR IR^. We end the paper with concluding remarks 
and discussion in section 6. Mainly in order not to interrupt the fiow of the paper with too 
many mathematical details we have written Appendices A, C and D. In Appendix B we 
consider the case where the background geometry has a covariantly constant null Killing 
vector Ijl^ (these geometries do not belong to the class considered already in section 2). 



2. General results 

Let us consider the (i-dimensional spacetime described by the metric 

ds'^ = 2 A{u,v) dudv + g{u,v) hij{x) dx^dx-' , (2-1) 

with = 1, 2, . . . , (i — 2). Let us also assume that there exist some matter fields with 
non-vanishing components of the energy momentum tensor given by 

T = 2 Tuviu^v.x) dudv + Tuuiu,v,x) du^ + Ty^{u,v,x) dv"^ -\-Tij{u,v,x) dx^dx^ . (2.2) 
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Notice that this form of energy momentum tensor is consistent with the Ricci tensor for 
(p|) as given by (A. 2) for / = 0. 

We consider a massless particle located at u = and moving with the speed of light 
in the f-direction and we want to find out what its effect is on the geometry described by 
( prrp . Similarly to [H our ansatz will be that for u < the spacetime is described by ( pTT| ) 
and for u > by ( |2.1| ) but with v shifted as v v + /(x), where f{x) is a function to be 
determined. Therefore the resulting spacetime metric and energy momentum tensor are 

(is^ = 2 A{u, V + e/) du{dv + Qf,^dx') + g{u, v + 6/) h,j{x) dx'dx^ , (2.3) 

where O = 0(w) is the Heaviside's step function audi 

T =2 Tuv{u, V + O/, x) du{dv + Of idx"^) + T^uiu, v + G/, x)dv? 

(2.4) 

+ Tyy{u, V + O/, x) {dv + Qf^idx"^)"^ + Tij{u, v + G/, x) dx'^dx^ . 
In order to compute various tensors it would be easier to transform to the new coordinates 



u = u, V = V + f{x)Q{u) , x' = x\ (2.5) 



in which the metric ( p.3|) and the energy momentum tensor take the form 

ds'^ = 2 A dudv + F dv? + g hij dx^dx^ 
F = F{u,v,x) = -2 A f 6 , 

and 



(2.6) 



T = 2 (f^^ - fyJS) dudv + (T^^ + fyypS^ - 2%yf5) dv" + f^^ dv"" + T,, dx'dx^ , (2.7) 

where the hats indicate that the corresponding quantities are evaluated at -u, v, and 
where 5 = 5{v) is the 5-function. 



For spherical sourceless shock waves in Minkowski space obtained with a different ansatz 



than (U) see Ij. 
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The various metric and field components must satisfy Einstein's equations in tlie 
presence of matter, which in d-dimensions read (throughout this paper c — h — G — 1) 

- ^a^uR - -Stt T^, R^, = -8% (T^, - g^, Tx^) = -Sir T^, . (2.8) 



Obviously T = T^j^dx^dx^ is also of the form ( p.2[ ). We would like to find solutions of 
( p.8|) with metric tensor of the form ( p.6|) and an energy momentum tensor given by ( |2.7| ) 
plus the contribution of the energy momentum tensor for a massless particle located at 
the origin of the transverse x-space and at u = and moving with the speed of light in 
the f-direction 

TP = du" = f?^ du" = -4pA^ 5(^^-2) (x) 6{u) du^ , (2.9) 

where p is the momentum of the particle. All the relevant tensor components appearing in 
( p.8|) are given in Appendix A. To simplify the notation we will also drop the hats over the 
symbols keeping in mind however the transformation ( p.5| ). Assuming that the parts of 
the equations ( |2.8|) that do not involve the function / are satisfied, one finds by examining 
the linear in f5 terms that at m = the additional conditions 

Q ^,11 T^-Lj 

d-1 n (2-10) 

A.,/ - -i^'-^f = P 9 A 5^'^-^\x) ^ 

must also be satisfied, where the Laplacian is defined as A/,,.^. = 1 / \/hdi\/hh'^^ dj . In order 
to cast the differential equation in (|2.1C1| ) into the given form we used the fact that at u = 



+ Sn T,, . (2.11) 



A 2 g 

This equality follows from the {fx, f) = (w, v) components of ( pT8| ) computed at w = and 
for / = 0. Notice also that the differential equation in ( p.lO|) does not explicitly depend 
on the components of the energy momentum tensor of the matter fields. Its dependence 
on these fields is only implicit through the functions A, g that are determined from the 
/-independent part of Einstein's equations. 



Next we examine the type terms. This is important because such terms should 
also vanish (in a distribution sense), otherwise our considerations are perturbative in pow- 
ers of /. Using the first line in ( p.lO|) it is easy to see that the coefficient of f'^d'^ in the 
Run component of the Ricci tensor in (A. 2) has terms of order 0{u) and O(w^) (in all of 
our examples such terms are of order 0{u'^) since there is functional dependence only on 
the product uv). Remembering that we are really considering all the quantities involving 
5-functions as distributions to be integrated over smooth functions we find that in this case 
such an integral vanishes. Moreover because T„„ = at -u = we have that Tyy = 0{u) 
(at least). Therefore we can take the terms in ( p.8| ) involving to zero. A last crucial 
remark is that because we have assumed that / is a function of the x*'s only, the potential 
V dependence in ( |2.1C1| ) should drop out for a consistent solution to exist. Mathematically 
that implies that the coefficient of the order 0{u) term in the expansion of g{u, v) in pow- 
ers of u should be a linear function of v. In fact this is the case in all the examples we 
explicitly work out. 

For the 4-dimensional case in the absence of matter fields the condition (|2.10|) was 
given in 0]. It is remarkable that (|2.10|) has essentially the same form in the vacuum and 
in the presence of matter except that in the latter case the additional condition T„„ = 
at « = should be imposed as well. Let us also emphasize that the differential equation 
in (|2.10|) is nothing but the Green function equation (of course with S^'^~'^\x) replaced by 
[x — x')). Thus instead of a point massless particle we could easily find the result if 
we had an extended source with density p(x), i.e. f{x) = J p{x)f{x, x')dx' . In this way we 
can consider spherical and planar shells of matter as in |jT^ or charged particles, cosmic 
strings, monopoles etc. (for the flat space case see [Q). 

The following remark is now in order. According to the work of in the 4-dimensional 
case in the equation for f[x) in ( |2.1C1| ) half (in our normalization) the 2-dimensional cur- 
vature R^"^^ constructed out of the metric hij should be in place of the factor In our 
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case using (|2.8| ) for (/U, z^) = the fact that g^y = at u = and (A. 2) we compute 



that at u = 

d-^ ^ ^ 4 f .^J (2.12) 

2 A 2 ^ ' 

Thus we obtain the result of [Q, and in fact generahzed in any number of spacetime 

dimensions, if Tijh^^ = or, after using (|2.8| ) ( pT2D , if T^v = 0. The stronger conditions 

Tij = and = 0, which imply Tij = and in our case T^v = were imposed in 0]. 

However, in many interesting cases, such as the Reissner-Nordstrom charged black hole, 

the component T^v is non-vanishing as we shall see. 

Although the explicit form of the metric, once the function f{x) is obtained by solving 

( p.lOp , is given by ( p.6| ) it would be helpful to examine the geodesic equations in order to 

obtain a clearer understanding of how the original geometry ( |2.1| ) is affected by the presence 

of a massless particle moving in the f -direction at w = 0. It is shown in Appendix C that 

as the geodesic trajectory crosses the null surface u = there is a shift in its t^-component 

given by 

Av = fix) , (2.13) 
as well as a refraction effect in the transverse x-plane expressed by the 'refraction function' 

Ri^) = ^\ n n+ = -| n/>'^ (2-14) 

that essentially measures how much the angle that the trajectory forms with the u = 
surface changes as we cross this surface. In other words a generic trajectory when 
it crosses u = suffers the discontinuity ( p. 13 ) in its f -component, with the u and 



components being continuous at that surface, and moreover its x^-components change 
direction along u according to ( p.l4|) . Although in both phenomena /(x) determines the 



functional dependence there is however a qualitative difference. In the case of ( |2.13| ) the 
discontinuity equals f{x) whereas in ( p.l4| ) the directional derivatives of f{x) play the 
important role. Therefore there might be points x* where there is a discontinuity in the 
trajectory but no refraction effect and vice versa. In fact this will be the case in some of 
our examples considered in the following sections. 



3. The (i-dimensional Reissner-Nordstrom charged black hole 

With this section we shall start applying the general formalism that we developed in 
the previous one by considering the case of a massless particle moving in the outer horizon 
of a Reissner-Nordstrom charged black hole in d dimensions. Since the latter is a spherically 
symmetric, asymptotically flat solution of (|2.8|) , our results will be a natural extension of 



those obtained by Dray and 't Hooft for the case of the 4-dimensional Schwarzschild black 
hole. The metric for the Reissner-Nordstrom charged black hole solution in d dimensions 



IS ^ 

ds"^ = -A(r) dt^ + A-^(r) dr^ + r^dflf^_2^ 

2C^ (3-1) 

^''^^ ~ j,d-3 ^2{d-3) ' 

where C, E are constants depending on the dimension of spacetime d and related to the 
mass and the charge of the black hole (see [|T2| ). The energy momentum tensor of the 
electromagnetic held is given by 



T,. = ^^^(^(^MAi^/ - \9,.F') , F,. = ^l(d-2)(d-3) , (3.2) 

where /i, = t, r, z and Ftr denotes the only no n- vanishing component of the electromag- 
netic field 2-form. It is possible to bring ( |3.1j ) into the form ( |2.1[ ) by means of the following 
transformation 

u = e'/'^F{r), ^ = e-*/«F(r), a = ^{l ~ { — (3.3) 

a — 3 r_|_ 

where r^~'^ = C ± VC*^ — -E'^, with r_|_ > r_, are the outer and inner horizons where the 
function A(r) in (|3.1| ) vanishes.i It turns out that 

F{r) = exp (— j drX~^{r)) = ^(r - r+)^ + . . . (3.4) 



Then since the only non-vanishing components of the electromagnetic field would be Fu 



the corresponding energy momentum tensor (3^) would be of the form ( |2.2| ) with Tuu = T^v = 0. 



Nevertheless, its precise form will not be needed for our purposes as we have already explained. 



where is a constant and the dots stand for higher order terms in the (r — r_|_) expansion.! 
We also find that 



1 2 f 

f ) = -Q!^A(r) exp ( / (irA~^(r)) , g 

Z ex J 



(^x,t;)=r2, (3.6) 



and that at w = the conditions g^^ = A^y = of ( |2.1C1|) are indeed satisfied. Using the 
values of the following functions at w = 

A^^n-\l-{-Y-')-\ g = r\, ^,.„ = 2r+^-^ (3.7) 
a — 3 ^ r_|_ ^ 



we find from (|2.10| ) that the equation the shift function satisfies is 

A(rf_2)/ - a(d,r+,r_) f = 2n b{d,r+,r.) d^''-^\x) , (3.8) 

where A(£;_2) is the Laplacian on the unit {d — 2)-sphere and a{d, r_), b{d, r_|_, r_) are 
constants defined as 

a{d,r+,r_)^hd-2){d-3){l-{-r-') 

^ _ (3.9) 

6(d,r+,r_) = 32prii(d-3)-^K-2(l-( — )'^-3) ^ . 

r_|_ 

Physically the most interesting is the 4-dimensional case. Then (|3.8| ) takes the form of 
(D.l) with 

J. ± i-a^ — 1-_ 



c=l (0<c<l), /c = 32pr*(r+ -r_) "+6 . (3.10) 

r_|_ 

Obviously if the black hole is uncharged (Schwarzschild), i.e. r_ = 0, we recover the result 
of Dray and 't Hooft . The solution of the equation is given in (D.2).i For i < c < 1 



^ The explicit expression for k is quite straightforward to write down in any number of 
dimensions. However the result is complicated and not very enlightening. For d = 4 the expression 
is simpler 



„2 



— r_i_ — r 

2 



K = (r+ - r_) e . (3.5) 

^ The solution of ( |3.8| ) for the higher than d = 4 dimensional cases can easily be given in 
terms of generalized spherical functions (see eqn. (D.12) in Appendix D). 
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(and in fact for all c > ^) an integral representation (proportional to a hypergeometric 
series) of the solution can be found using 

/•oo I A- - 

I '^-"'"*"'=-(v^^») = ^(J^. (3.11) 

and the generating function for Legendre polynomials 

oo 

^t^Pz(cos^) = (1 -2cos^ t + ^2)-^/^ . (3.12) 

1=0 

We find (consult 0) 



k f°° 

f{e-c) = --j=j^ dscos{^c-l/A s) 



1 



cosh s — COS 9 



F{- - i^^^ljl, - + i^^^ljl, 1, cos^ 



(3.13) 



2cosh(v/c-l/4 tt) '2 ^ ' ' 2 ^ ' ' ' 2'' 

For < c < I we should replace \Jc — 1/4 by i ^Jl/A — c with parallel replacement of the 
corresponding trigonometric functions by hyperbolic ones and vice versa. Notice that the 
solution blows up at the point of the unit 2-sphere where the particle was placed, i.e. at 
^ = 0. Moreover it is everywhere negative and for fixed c it is a monotonically increasing 
function of G [0, tt] approaching a constant at ^ = tt. For fixed 9 it also monotonically 
increases as a function of c G (0, 1]. The 'refraction function' is 

mc) = {^)^^,def{9-c) . (3.14) 

As a function of 9 it monotonically decreases from plus infinity at ^ = to zero at ^ = tt. 

It is useful to examine what happens in the region close to the two poles of the sphere 
where the behavior is extreme.! Close to the northern pole at ^ = we have 

-kKo(JE9), as^«l 
f{9-c)-{ ^ (3.15) 

k ln(v/c 9) , as 9 «< 1 , 



^ In the following In{x) and Kn{x) will denote, as usual, the nth-order modified Bessel func- 
tions. The relations: K'q = —K\, I'q = I\ and the leading order behavior as x << 1: /q — 1, 
/i ~ |, Kq ~ — Inx, Ki ~ ^ will also prove useful. 
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whereas close to the southern pole at 6* = tt 



=^ Io{^/c (tT - 6')) , as TT - 6' << 1 



2 cosh(-v/c— 4-7r) 

f{0;c)^{ ; (3.16) 
^^"^ , ^ , as TT - 6* «< 1 . 

2 cosh( y' c— ^71") 

For the 'refraction function' the corresponding results are 

( k' y/c Ki{y/d e) , ase«l 

R{e;c)c:^{ (3.17) 

[ |- , as 6' <<< 1 , 

and 

2 coshi \ c— 4-7r) 

i?(^;c)^<' (3.18) 



■KC 



=^(7r-6'), as7r-6'<<<l 



4 cosh(-^ c— ^7r) 

where /c' = (A/(7)^=o^ is a positive constant. As it was expected as we approach 6* = 
both / and R behave like the corresponding functions in the flat space case 00 (see 
also (B.8)) and both approach infinity. At the southern pole both functions reach their 
minimum magnitudes and in fact the refraction phenomenon disappears even though a 
particle trajectory is still discontinuous since /(tt; c) 7^ 0. 

We end this section with a comment on the case of the extremally charged Reissner- 
Nordstrom black hole. It is obvious either from the series representation (D.2) for f{9;c), 
or from the integral one ( |3.13| ) that for c = 0, or equivalently in the extremal limit r_ = r_|_, 
the solution is not well defined. In fact the choice we have made for the constant a in 
( p.3| ) and everything that follows are not valid if r_ = It can easily be shown that 
there is no solution to our problem, or in other words no single particle can move with the 
speed of light in the outer horizon of an extremally charged Reissner-Nordstrom black hole. 
Physically this should have been expected because extremallity is the situation where the 
two horizons coincide just before a naked singularity appears. From a more mathematical 
point of view this happens because the function F^(r) is not analytic at r = r+ = r_, as 
it is in the non-extremal case (see (|3.4| )). 
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4. Constant curvature spaces 

In this section we consider the De-Sitter and Anti-De-Sitter spaces where there is 
a constant curvature R = 2d/{d—2)A corresponding to a non-vanishing cosmological 
constant A. Even though a cosmological term was not considered in ( p.8| ) its effect can be 
imitated by an energy momentum tensor in ( p.8|) that is proportional to g^,^, i.e. 8nT^,^ = 
—2/{d — 2)Kg^jj. We also analyze the Schwarzschild-De-Sitter black hole case. 

4.1. The 4- dimensional cases 

The metric for the 4-dimensional De-Sitter space is 

2 2 

dsls = -(1 - ^) dt^ + (1 - ^)~^dr^ + r\d9^ + sin^ 9 dcj)^) , (4.1) 



whereas for the corresponding Anti-De- Sitter one is 

2 2 

dslos = - 1) dt^ + ly^dr^ + r^{dx^ + sinh' x d(t>^) . (4.2) 

In both cases the constant a is related to the non-vanishing cosmological constant!^ as 
A = ±3/a2. To bring (|T|), (gj) into the form (|T|) we make the same coordinate 



transformation as in ( |3.3| ) but with a = a. Then it turns out that 

F(r) = (±^)^, Aiu,v) = l{r + af , g{u,v)=r^. (4.3) 
r + a 2 



The conditions of (|2.10|) are indeed satisfied at w = (r = a) and the differential equation 



/ satisfies is of the form (D.l) with 

c = ^ 2 , k = 32pa^ . (4.4) 
Then according to (D.5) for / = 1 the solution for the De-Sitter case is 

fosiO) = 32pa^ (1 - ^ cos^ ln(cot2 ^)) e(| - 6) , (4.5) 



^ Throughout this subsection upper (lower) signs in some quantities correspond to the De- 
Sitter (Anti-De-Sitter) case. 
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where contributions proportional to Y{^{9,(f)), m = —1,0,1 that solve the homogeneous 
equation have been omitted. As it also pointed out in the Appendix D the use of the 
0-function that essentially restricts the solution to the upper hemisphere is necessary in 
order to have a solution that does not blow up at the southern pole at ^ = tt, which 
would be unphysical since no particle is placed there. Comparing with the black hole case 
we considered in the previous section we find a notable difference. The solution goes to 
minus infinity at ^ = and then monotonically increases until it reaches the value (32pa^) 
at ^ = ^. Therefore there is an angle (^o — 33. 52^^) where it becomes zero and then it 
passes from negative to positive values. Thus contrary to intuition the magnitude of the 
shift has its minimum not at ^ = tt (as it was the case in the previous section) but at an 
intermediate angle. The corresponding 'refraction function' is 

RDsie) = 6W + ^ sin^ ln(cot2 ^)) G(| - 6) - Mpa^ 5(9 - |) . (4.6) 

The first term is a monotonically decreasing function of 9 and it varies from plus infinity to 
zero as we go from the northern pole to the equator. However exactly there the second term 
gives an infinite contribution. This is in effect a consequence of restricting the solution to 
the upper hemisphere and, in some sense, can be thought of as a source term that displaces 
the minimum magnitude of fosid) from 9 = ^ to 9 = 9o. Notice also that at 6* = 6*0 there 
is no discontinuity for the geodesic trajectories but there is still a refraction effect. 

A natural question to ask is whether or not there is a distribution of massless particles, 
instead of just a single one, for which the use of a 0-function in (4^) is unnecessary. Ob- 



viously there are many such distributions (after all (|4.5|) is nothing but a Green function) . 
The simplest one is to consider two particles of the same energy p one in each pole of the 
2-sphere. Then the solution becomes exactly i\i.5\) with no 0-function of course, and in 
fact it can be found by infinitely boosting a Schwarzschild-de-Sitter black hole (see (4.9) 



below) in the zero mass limit \10 
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For the Anti-De- Sitter case (D.8) for / = 1 gives 

fADsix) = S2pa^ il coshx ln(coth2 |) - l) ■ (4.7) 
Notice that in the Anti-De-Sitter case we did not have to 'cut' the solution by itroducing 



a G-function because (|4.7|) vanishes by itself for large x- Of course this is due to the 
fact that in contrast to the De-Sitter case the 2-dimensional space the shift function / 
takes values on is a non-compact hyperboloid and not the 2-sphere. The corresponding 
'refraction function' is 

RADsix) = (-|^ + I sinhx ln(coth2 |)) . (4.8) 

As a function /ad six) {Rads{x)) monotonically decreases (increases) from plus (minus) 
infinity to zero as we get away from the pole of the hyperboloid at x = 0. This is similar 
(but not exactly the same) to the behavior of the corresponding functions in the previous 
section. 

Let us now return to the more general case of the Schwarzschild-De-Sitter black hole 
for which also there a positive cosmo logical constant. The metric is 

dsljys = -(l- — -^)dt' + (l- — - ^y'dr^ + r\de^ + sin2 6 d<P^) . (4.9) 

Since it is a straightforward calculation we decided not to give algebraic details concerning 
the derivation of the various results. The differential equation to be satisfied by the shift 
function f{9) is again of the type (D.l) and it turns out that depending on the ratio a/M 
there are two branches of solutions for the constants c and k. We will only consider the 
branch in which the null surface where we will place the massless particle corresponds to 
a positive value of r (for positive M). In this branch 

in r- (f 09 , — M 

c = 2 sm — ( V 3 cos sm — ) , cos (p = v 27 — 

3 3 3 a 

^>V27 [-|,o)u(o,^] , ce [-2,o)u(o,i] , 
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and the constant k is always positive with precise value which is not of any particular 
interest. The null surface u = where the massless particle is placed corresponds to 
r = rs = 2^cos^^^. For positive (negative) values of (p the radial coordinate r lies 
between < r < ri (ri < r < rs), where ri = 2^ cos { ri, and r2 = — 2^ cos ^ 
are the three surfaces where gtt vanishes). The boundary cases 99 = ^ corresponding 
to the De-Sitter space (c = —2) and the Schwarzschild black hole with zero cosmological 
constant (c = 1) have been already discussed. The case = = v^27, c = 0) 
is excluded because the situation is similar to the extremal Reissner-Nordstrom charged 
black hole. Except for the case c = —2 the solution is given by (D.2). For the range 
< c < 1 the behavior of the solution is exactly the same as in the case of the Reissner- 
Nordstrom charged black hole of the previous section. For the range —2 < c < 0, similarly 
to the case of ( |3.13|) , we find the following integral representation 

f(e;c) = -^-k r dscosh{,/T/I^cs) ( . - e-/^) . (4.11) 

Jq Vcoshs — cost' 

Comparing with (|3.13| ) we see that there is a second term in the integrand that essentially 
regulates the divergent behavior of the first term for large values of s. We note that the 
solution again blows up at ^ = and as in the case of the De- Sitter space it is monotonically 
increasing as we move from the northern to the southern pole of the 2-sphere and changes 
from negative to positive values at an angle ^0 that depends on the value of c. For instance 
for c = —1 we have 9o — 61.59^. Thus similarly to the De-Sitter space case the magnitude 
of the shift function is reaching its minimum at a point which is not the furthest from the 
massless particle at ^ = 0. The important difference with the De-Sitter space case is that 
the solution is extended to the southern hemisphere as well, i.e. no use of G-functions 
is required. Let us also note that the 'refraction function' is a monotonically decreasing 
function of 9 and that the functional dependence in the region around the northern and 
the southern poles is again given by ( ^.15D -( |3.18[ ) (with a different value for k and where 
the appropriate analytic continuations in the hyperbolic functions should be performed). 

Let us also mention that we refrained from presenting here any results concerning 
analogs of Reissner-Nordstrom charged black holes (see [ pO| ) with a cosmological constant 
because that would have complicated things rather unnecessarily. 
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4-2. The 3- dimensional black hole 

It is interesting to consider the Anti-de-Sitter space in d = 3. Then the metric is 

2 2 

ds^ = -(^ - 1) dt^ + - ly^dr^ + r^dx^ , (4.12) 

where x is a non-compact coordinate. Again the constant a is related to the cosmological 
constant which is negative, i.e. A = — 1/a^. The change of variables that brings ( |4.12| ) 



into the form (|2.1| ) is similar to (^4.3| ) with the lower sign in the expression for F{r). Then 
the 1-dimensional equation that / satisfies using ( |2.10|) for d = 3 turns out to be 

^ - / = 64npa^ d{x) ■ (4.13) 

Its solution that has the correct asymptotic behavior reads (we ignore the obvious solution 
of the homogeneous equation) 

f{x) = ~32npa^ e"!^! . (4.14) 

A more physical situation arises when the non-compact variable x is made a compact 
one by letting x — ^ and identifying (p = (p + T. The most physical choice is of course 



T = 2tt. In this case we have a black hole in (i = 3 as it was shown in In fact this 
solution can also be thought of as a solution to string theory with non-trivial antisymmetric 
tensor and constant dilaton fields [^. The corresponding exact conformal field theory is 
obtained if one quotients the 5'L(2,IR) WZW model by a discrete subgroup. Because of 
the above discrete identification the 5-function that appears in ( [1.13|) should be replaced 
t'y J2^=-oo ^i'P ^ nT). Then instead of ( |4.14[ ) the appropriate periodic solution satisfying 
/((/> + T) =/((/>) is 

oo 

/((/)) = -32npa'^ ^ g"!'^-"^! , |(^| < oo . (4.15) 

n= — oo 

After some algebra the infinite sum is computed to give 

f{<P) = -GATvpa^ (-^ cosh<^-f ie-l^l) , \<P\<^- (4-16) 
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Notice that the range of is restricted in the fundamental domain and that for T — oo 
we recover the expression ( [4.14|) as we should. On physical grounds we expect that the 
absolute value of the shift is always larger that the one predicted in the non-compact case, 
for any (f). This should be the case because the identification (p = (p + T 'creates' many 
different sources of infinite curvature (5-functions) in the entire real line. Using ( [4.16|) (|4.13|) 



one easily verifies that statement. In fact |/(^)| reaches its maximum (minimum) value at 
(p = {\(p\ = T/2). Finally, the 'refraction function' is 

R{<P) = -1287rpa4 sinhc^ - hign{cf>)e-\^\) , | ■ (4-17) 



5. String inspired gravitational solutions 

In low energy heterotic string theory we may consider the set of background fields 
that comprises a metric, an antisymmetric tensor, a dilaton and an electromagnetic field 
coupled in a 2-dimensional a-model action. Then the requirement of conformal invariance 
generates constraints these fields have to obey, the so called beta-function equations. For 
our purposes it is convenient to present these equations in the Einstein frame where the 
metric gj^^^ is related to the cr-model one by 

2 

(7^L =exp(^— ^$) (7^, . (5.1) 

In the rest of this section we drop the extra index having in mind that all quantities are 
evaluated in the Einstein frame. To lowest order in the string coupling constant a', the 
beta-function equations (see for instance ||2^) arell 



R>^^ = ^{D^<5D,<5 - g^.D^<5) - -^{F^^F,^ - ^9^..F^) exp(^$) 
2$ 

^A(exp(^)F/)=0. 



(5.2) 



^ We will not mention the dilaton beta function since it is not independent from the other 
ones. Also since in the examples we will shortly give the antisymmetric tensor is locally trivial, 
i.e. Hf^up = 0, we restrict to only such cases although it is straightforward to workout out the 
details in the more general case as well. 
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The first of tlie above equations is of the form of Einstein's equations (2^) with non-trivial 
matter energy momentum tensor that is being created by $ and F^i, and given by the 
right hand side of this equation. In the following we consider metric tensors of the form 



given by and electromagnetic field strength and dilaton field of the form 

F = 2 Fuviu, v) du Adv + Fij{x) dx" A dx^ 

(5.3) 

$ = v) . 

Strictly speaking with the above form for the electromagnetic field -F^y the antisymmetric 
tensor field strength if^j^p cannot be set to zero since its source term dH = —F A F does 
not vanish. However, one may think of having two independent electromagnetic fields 
F", a = 1, 2 given by 

F^ = 2 Fl^ (m, v) duAdv , = F^. (x) dx' A dx^ . (5.4) 

In that case dH = —F"- A F" = 0. The two formulations are equivalent since in both cases 
the electromagnetic energy momentum tensors are the same leading to the same equations 
of motion given by ( ^.2|) . 

As in section 1 we 'perturb' a given background solution by adding a massless particle 
moving along the f-direction at w = 0. Since the first equation in ( |5.2|) is the same as 
Einstein's equation ( p78|) the conditions on the functions g(u,v), A{u,v) and f{x) are 
precisely given by ( |2.1C1| ), whereas T„„ = at m = translates into a condition on the 
dilaton (again at m = 0) 

= . (5.5) 

Examination of the second equation in ( p.2|) gives the condition {Fuv),v = at m = which 
however is obeyed thanks to the equations of motion and the forementioned conditions on 
A{u, v) and g(u, v). 
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5.1. An electrically and magnetically charged black hole 

As a first application of tfie above let us consider a solution to low energy heterotic 



string theory that represents a black hole with both electric and magnetic charges [jl^ [|T5 



dsl = -Xi{r) dt^ + X^\r)dr^ + (r^ - r^) {d9^ + sin^ 9 d(f>^) 

^^^"(^TT^) (5-6) 
Qe 



where 



,^ dt Adr + Q m sin 9 d9 A dcj) , 
- roy 

_ {r-r+){r-r-) _ Qm - Ql 



(5.7) 

r±=M±^M^ + rl-Ql-Ql, . 
As in the case of r_(_ and r_ are the outer and inner horizons respectively. The 

constant tq is the so called dilaton charge that essentially measures the difference between 
electric and magnetic charges. One easily finds that a transformation of the type ( |3.3| ), 
but with a = 2(r^ ~ ^o)/ (^+ " ^-): takes the metric into the form (|2.1| ) with the relevant 
functions given by 

2 2 

F(r)=e'^/-(r-r_)"^<^T^(r-r+)^ , A{u,v) = ^Xi{r)/F{rf 
g{u,v) =r'^ -rl . 

Obviously after this change of coordinates the expressions for the dilaton and the elec- 
tromagnetic field strength assume the form ( |5.3|) . One easily finds that (|2.10|)(|5.5| ) are 
satisfied at u = (or r = r_|_) and that the equation for the shift function / is given by 
(D.l) with 

2 _ 2 

^ ^_ r_|_ — T 

c = r JX~''~ , k = ?,2p{r\-rlf{r+-r_)~^+^e~''^^+^- (5-9) 

Similarly one easily sees that because of the fact that tq < r_ < r_(_ we have < c < 1 
and therefore the solution to / is given also by the general expression (D.2) or (|3.13| ). 

19 



5.2. The coset 31(2,11) /JR 011^ 

Finally let us consider the case of the 4-dimensional model that is the tensor product 
of the 2-dimensional coset 5'L(2, lR)_jt/IR |l2^] with 2 additional non-compact dimensions. 
The metric in the Einstein frame (c.f. ( |5.1| )) and the dilaton are 

ds'^ = — 2 e dudv + (1 — uv){dx^ + dy"^) 

(5.10) 

$ = ln(l - uv) , 

with e = sign(/c), where {—k) is the central extension of the 5'L(2,IR)_fc current algebra. 
For e = 1 the causal structure of the spacetime is that of a black hole with a singularity 
at future times t = u -\- v, whereas for e = — 1 it has the cosmological interpretation of an 
expanding Universe with no singularity at future times t — u — v (see [^). The metric is 
already of the form (|2.1|) with 



A{u, f ) = — e , (/(m, v) = 1 — uv . (5-11) 

It is easy to see that ( |5.5| ) is satisfied and that the differential equation for the shift function 
is 

' + - e)/(p) = -16ep -5{p) , (5.12) 



dp^ p dp p 
where p^ = -\- y'^. 

Let us first consider the black hole case (e = 1). Then ( |5.12| ) is nothing but the 
modified Bessel equation of order zero with the additional 5-function source term in the 
right hand side. Its solution with the correct asymptotic behavior is 

/+(p) = 16pKo(p) 

(5.13) 

i?+(p) = 16pKi(p) . 

In order to obtain a 5-function behavior in the right hand side of ( p.l2| ) one should carefully 
treat the derivative part of the left hand side using the fact that -K'o(p) — ~ ln(p) for p << 1. 
Also if one wishes to compactify one or both of the a;, y coordinates one should replace the 
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solution ( |5.13|) by a periodic one that will necessarily involve infinite series as it was done 
in the case of ( |4.14| ) and ( ^4.15| ) . 



In the cosmological case (e = —1) (|5.12| ) becomes the usual Bessel equation with 



solution respecting the singular behavior of the source term in the right hand side given 
by 

=87rpiVo(p) 

(5.14) 

R.{p) = -Sivp Ni{p) , 
where Nn{p) denotes Bessel functions of the second kind.i Obviously ( |5.13| ) ( pTl^ ) are 
similar to (|3.15|)(|3.17| ) since in both cases the metric in the transverse x-plane is the flat 
one (notice: No{p) ^ ^ln(p) and Ni{p) ~ for p << 1). However, whereas in the 

black hole case the shift and the 'refraction function' never change sign, in the cosmological 
case they do so inflnite number of times until the die off completely at p = cxd. This type of 
behavior most likely will show up in other cosmological models as well. The difference in 
behavior is more dramatic far away from the particle's position at p = since for p >> 1 
^o(p) - \Jf^^~^ whereas A^o(p) - y^sin(p - |). 

6. Discussion and concluding remarks 

In this paper we have found the necessary and sufficient conditions for being able to 
introduce a gravitational shock wave in a quite general class of background solutions in 
Eistein's general relativity coupled to non-trivial matter (with or without a cosmological 
constant) and string theory. We have also applied the general formalism to many cases that 
are of particular interest such as charged black holes and constant curvature spacetimes. 
The are also a few points we would like to further discuss. 



^ A more systematic way to obtain the solution of ( ^.12| ) is to expand as f{p) = 
dkkA{k)Jo{kp) and use the fact that = dkkJo{kp). Then easily we obtain A{k) = 
. Performing the integral we get for the shift function the expression in ( |5l3|) or (|5l4|) 
depending on the value of e. 
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1. It is quite interesting to notice that the addition of the massless particle in a 
background geometry creates in effect a perturbation that can be described, in the context 
of string theory in curved spacetime, in terms of a massless vertex operator. To be precise 
consider Einstein's equations in the vacuum or in a bosonic string theory language a 2- 
dimensional cr-model with zero antisymmetric tensor and dilaton fields. Let's take 

V ^ F^,{X)dX^dX'' , (6.1) 

as a candidate for a massless vertex operator.0 If we want to perturb the cr-model by this 
operator we should require that in every order in perturbation theory in the string coupling 
constant a', its anomalous dimension is 2 so that the perturbation stays marginal. That 
gives a set of consistency conditions F^^, is constrained to satisfy. To leading order in the 
string coupling they read 

- D^Ff,, - D^D^Fx^ + R/%Fp^ + D^^D^F^^) = 

(6.2) 

D^Fx^ - DPD'^Fp^ + RP'^Fp^ = . 
It is a straightforward (though a bit lengthly) calculation to verify that 

Fuu = -2A{u,v)f{x)5{u) , (6.3) 



as it is read from ( |2.6| ), is indeed a massless vertex operator corresponding to the back- 
ground metric ( |2.1[ ) provided that, as in ( |2.1UD , at w = we have g^y = and A^y = 



and that f{x) satisfies the homogeneous differential equation in (|2.1C1| ). The vanishing of 
the right hand side occurs because we have not included the energy momentum tensor 
corresponding to this vertex ( ^.91 ). This is not usually done if the background is flat be- 
cause then the resulting homogeneous equation has a solution (with a different anzatz for 
a vertex operator (see [|2^)), which is not however always true, for instance, as we have 
seen in the cases of black holes. 



In general j25[ j one may add the term a'R'^^^F{X), where i?'-'^-' is the curvature of the 2- 
dimensional worldsheet and F{X) is a tachyon-like operator. However this is not necessary here. 

22 



2. What is the meaning of the diverging solution (see ( |3.13|) ) as we go to the extremal 
limit of a Reissner-Nordstrom charged black hole and moreover the no solution result we 
have found working exactly at the extremal limit? In general the existence of a solution 
means that the original background responds to the 'perturbation' created by the massless 
particle. A particle geodesic, in a charged close to the extreme limit black hole, gets 
shifted and refracted by large amounts due to the large gravitational field created by the 
massless particle. Exactly at the extremal limit there is no solution with a single particle 
but there is one if in addition to the positive energy massless particle at the northern pole 
we place another one with negative energy of the same magnitude at the southern pole 
(both at r = r_|_ = r_). Then the solution with the correct singular behavior at ^ = 0, tt 
is f{6) = —kQo{x) = — ^In^cot^ |). In other words we might think that as r_|_ r_ 
and the gravitational field becomes infinite it is necessary to take the antiparticle into 
consideration as well. The fact that there is a solution in the case of an extremal Reissner- 
Nordstrom black hole only if we use a particle and its antiparticle is somewhat analogous 
to the 'semiclassical' explanation of Kleins's paradox using the concept of a 'sea' filled with 
negative energy particles. 

As a final comment on the Reissner-Nordstrom charged black hole let us mention what 
the result would be if we had placed the massless particle in the inner horizon instead of 
the outer one. The differential equation for the shift function is again of the form (D.l) 
with 

2 2 
—r 

c=l-— (-cx)<c<0), k = -32prt{r+ - r_) . (6.4) 

r_ 

Thus c, k are basically given by ( |3.10| ) with r+ and r_ interchanged. The solution for the 
shift function is given by (D.2) which one can show is equivalent to 

n — l I 4- - 1 / /2 

Ue) ^-kJ2 ^^(^o^^) - M cosh(VT74^ .) { . 

^^l[l + L) + c Jq V cosh s — cosy 

n-1 

- e"(^+i)' Pi{cose)) , -n(n + 1) < c < -n(n - 1), n = 1, 2 . . . . 

1=0 

(6.5) 
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Using (|6.5| ) one finds that the shift function starts from plus infinity at ^ = becomes 
zero and changes sign as many times as the value of n and eventually assumes the finite 
constant value predicted by ( |3.16|) . Notice also that ( [4 .111 ) corresponds to the case n = 1 of 
the above formula. Of course for the cases where c = — r2,(?i + l), ti = 1, 2, . . . (D.5) or (D.6) 
gives the solution. The case c = — 2 (ti = 1) is exactly the same as in the De-Sitter space 
and the solution is unique. For higher values of n however, there are (n — 1) undetermined 
coefficients which can probably be fixed with some additional physical requirements for 
the shift or the 'refraction function'. 

3. Can we find the shock wave-type of disturbances created by two or more massless 
particles moving along parallel surfaces (different values of u) just by a coordinate shift? 
Obviously this requires A^y = g^y = at these values of u. This is not possible for black 
hole type solutions. It can however be done for backgrounds with a covariantly constant 
null Killing vector (see Appendix B) that include plane waves and the fiat space as special 
cases. Then at each null surface we have a shock wave with a shift function obtained by 
solving (B.6). 

Finally, it would be obviously very interesting to explicitly compute, using techniques 
similar to those in 000, particle scattering amplitudes in the shock wave curved geome- 
tries we have derived. In particular one would like to know how these amplitudes depend 
on the different behaviors of the shift and the 'refraction function' we have demonstrated 
in the present paper. 
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Note added 



Towards finishing typing this paper we become aware that the authors of ||2^ dealt 
with some of the issues of this paper. However, we disagree with many of their findings. 
The main points of disagreement are (equation labels starting with 'LS.' refer to |P7|| ): 

1) For the general (i-dimensional case we do not agree with the expressions for the Ricci 
tensor (LS.5-9). Our (A.2) contains some additional terms with explicit d dependence. 

2) We do not quite agree with the conditions (LS. 17-18) (we find instead ( p.lQ )) mainly 
because of the explicit appearance of the cosmological constant in (LS.18) and the absence 
of a condition on T^^ in (LS.17). 

3) For the case of a Reissner-Nordstrom charged black hole (LS.35) agrees with (D.2) 
(with a = c). However, we find the different integral representation ( |3.13|) instead of 
(LS.36). 
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Appendix A. Components of useful tensors 



The corresponding to the metric (|2.6| ) non-vanishing Christoffel symbols are (for no- 
tational convenience we omit the hats) 

-pU _ I ■^,U -pu _ 9,V 7 

^uu- 2A^ A ' 2A 

F FF FA F F- 

~ 2A 2^2 ' ~ 2A ' ~ 2A 

™~ A ' ~ ^ 2A 2^2 ' ^ 

pi __}_ hik p -pi 9 ,u c-i -pi 9 ,v tri 

= -jh''\hik,j + hij^k - hjk,i) . 

Using the above expressions we find that the non- vanishing components of the Ricci tensor 
are (we substitute F = —2Afd (see (|2.6p)) 



p ~ ^ { 9,uA^u 9,uu A r . r d-2 g^y , 

+ (2 % - 2 + ^ (^,.A,„ + g,yA,^)) S f 

+ ^ ^ A A2 2 gA ^ 

P f-A,u-A,v -A,uv I ^ 2 g.u9,v d 2 g.uv\ 

, ('A,v A^yy d 2 g^yA^y^ 
^ " ^4 2 ^ 

P 2 /g^yA^y g^y 9,vv\ 



(A.2) 



In the 4-dimensional case ( |A.2| ) reduce to the Ricci tensor given in 0. In the case we are 
dealing with string theory as in section 5 it might be useful to compute the tensor Df^D^^ 
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with $ = f ). Its non- vanishing components read 



D D ^ = ^ 



A 



(A.3) 



Appendix B. Backgrounds with a covariantly constant null Killing vector 

Consider the metric with a covariantly constant null Killing vector 

ds^ = 2dudv + gij{u,x) dx^dx^ , 

where = l,2,...,(i— 2). The most general matter field energy momentum tensor 
consistent with the non-vanishing components of the Ricci tensor corresponding to ( P . 1| ) 
(see (B.5) below) has the form 

T = Tuu{uj x) du^ + 2 Tui{u, x) dudx^ + Tjj (u, x) dx'^dx^ . (B.2) 

The analog of ( pT6|) is 

ds^ — 2 dudv + F dv? + cjij dx^dx^ , 

(B-3) 

F = F(u,x) = -2 f S , 

whereas that of ( |2.7D is exactly the same tensor ( [B.2| ) since no t^-component of the energy 
momentum tensor appears in ( p.2| ). 

The non-zero components of the Christoffel symbols and the Ricci tensor are (once 
again we drop the hats and also we denote derivatives with respect to u with a dot) 



^uu r, 1 ^ ui r, ' ^ij 9ij 

^uu = -77 9' Fk , • = -g' gjk , = -g' (gi^j + gij,k - 9jk,i) , 
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and 

Run = {-^9ij9'^ + ^9'^9^^9ij9ki) + S Ag^J 

Ru^ = -D^^gk]j9'^ (B-5) 

Then the condition to be satisfied at u = is the difi'erential equation (again we drop 
the hats over the symbols) 

Ag^J = 32n p , A,,^. = ±d,g^^ ^d, . (B.6) 

V9 

Let us next consider string backgrounds with a covariantly constant null Killing vec- 
tor, i.e. with metric in the Einstein frame as in ([B.l|) , and in addition with non-trivial 
antisymmetric tensor field strength components Hijk = Hijk{u,x), Hiju = Hiju{u,x) as 
well as a dilaton field $ = x). It can easily be shown that the 'effective' matter energy 
momentum tensor due to these background fields is again of the form ( p.2|) . Thus, the 



equation for the shift function is again given by (|B.6|) , where the metric is the Einstein 



frame one. For completeness, the non- vanishing components of -D„-Di^$ are 



D^D,^ = - ]^9^'gH^,k (B.7) 

Let us also note that for the particular case of plane waves where all scalars and tensors 
depend only on the light cone coordinate u and not on the x*'s, equation ( |B.6| ) is the same 



/WH ,,,,,, (B.8) 



D/2 



as the one corresponding to the fiat space case and therefore its general solution is 0] 

16p In p , if d = 4 

327rp 1 

(d-4)nd_2 p'*- 

where in general = r{D'/2) 'area' of the unit sphere in D-dimensions and = 

gij{0)x^x^ . Since we can freely shift w by a constant the choice of u = should be made 
in such a way that the metric is non-degenerate at that point, i.e. m = is not a focusing 
point of null rays. 
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Appendix C. Geodesic equations 

For the metric ( |2.6| ) the geodesic equations obtained by varying v and are (in this 
Appendix the dots over the various symbols denote derivatives with respect to the afRne 
parameter r) 



I ^i"" ■ 2 9,v , ■ i ■ j , p -^,v r • 2 n 

U + U — — ^ hijX X-' + / W = 

A 



x' + T\x^x^ + ^ ux' + ^ vx' + - 5 fM' = 
whereas the one obtained from the variation of u is 



(C.l) 



I -^^v ■ 2 9,u 7 . i . j 

V + — — V — hiiX X-' 
A 2A ^ 

+ (/^ - 2/^ iiv - 2/, ux^ - ^-^f h^jx^x^) S (C.2) 

- f 6' + 2f 6^ ^ ii^ = . 

A 

Prom (|2.6| ) the energy corresponding to the geodesic is given by 

2Auv + gh.jx'x^ - 2f5Au^ = a , (C.3) 

where a = —1, 0, 1 depending on whether the geodesic is timelike, null or spacelike respec- 
tively. It is clear from ( |C.3| ) that the it^-term can be compensated by a discontinuity in v 
at w = 0. Taking v = vq + 0{u)Av, where vq is a solution of the geodesic equations for 
/ = 0, and integrating ( |C.3| ) over a small interval around w = we obtain immediately 
the relation (|2.13|) . Performing the same integration in the first of (|C.1J ) gives no new in- 



formation because A^y = at m = 0. However the same procedure in the second equation 
in ( |C.1|) gives 



A, 



U=0~ \U=0+ g 



Our next step is to express the variation with respect to the affine parameter r at m = 
as a variation with respect to u itself. This can be done by noticing that the first equation 
in ( |C.1| ) evaluated at w = is 
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and has as its solution 

Therefore ( |C.4|) takes the form of ( p.l4|) . It is important to emphasize that ( |C.6|) is reaUy 
a solution of (|C.5|) only at u = where A „ = 0. Finally, integrating ( |C.2| ) over a small 
interval around w = gives no additional conditions as it eventually reduces to ( |(J.5[ ). 



Appendix D. Solutions of the equation for the shift function f 

Let us consider the differential equation 

A(2)/ - cf = 2Txk 5{x - l)5{4>) 

Q2 (D.l) 

A(2) = -a;^)^^; + - — a; = cos 6', 
^ ' 1 — 

on the unit 2-sphere with metric ds^ = d6^ + sin^ 6 dcfP , where /c, c are real constants. 
Except for the 5-function in the right hand side this is nothing but the usual Legendre 
equation of order where v is a. solution of z^(z^ + 1) + c = 0. Therefore its solutions 
depend very much on the value of the constant c. 

We first consider the case of c 7^ — n(n + 1), n = 0, 1, . . .. Since the eigenvalues of the 
Laplacian operator are = —n{n + 1), n = 0, 1, . . . the modified Laplacian ( A(2) — c) has 
no zero modes and therefore it is invertible. The general solution of ( p.l| ) can be given in 

terms of Legendre polynomials as 

00 ^ _^ 1 

f{e; c) = -kJ2 Plicos 9) , c e IR - {-n{n + 1) , n = 0, 1, . . .} . (D.2) 

z=o ''^^ + -LJ +c 

Since in this case c does not coincide with an eigenvalue of the Laplacian there is no 
solution to the homogeneous equation in ( p.l| ). 

If c = + 1), I = 0, 1, . . ., then we should make use of the eigenfunctions of the 
Laplacian that are not regular at the poles at a; = ±1, namely Qi{x) that are defined as 

Qi{x) = }-Pi{x) In(l^) - -Pm-i{x)Pi-Ux) . (D.3) 

m—1 
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If the solution for / is proportional to such a function then a careful treatment of the 
derivative part in ( p.l|) will produce 5-functions at the poles of the sphere at x = ±1. 
However according to the right hand side of ( p.l|) there should not be any singularity at 
X = —1. Therefore we are forced to consider a solution of the form / ~ Qi{x) G(x — xq), 
where xq is to be determined by requiring that / satisfies ( p.l|) . After some algebraic 
manipulations and using the fact that -^0(x) = 5{x) we obtain the condition 



{1 - x^)^^^^ S{x - xo) = . (D.4) 



Therefore xq must be one of the points where the slope of Qi{x) vanishes. In fact there are 
/ such points {xi, X2, . . . x^}, symmetrically distributed around x = which always satisfies 
( p.4|) if / is odd. Thus the general solution is even by 



i-i 

2 

f{9) = -k Qi{x) (Aoe(x)+ ^[A^e(x-x^) + S^e(x + x^)]) , x = cos^ 

m=l 

l-l 

2 

Ao+Y,i^rn + B^) =1 , / = 1,3,... 

m=l 

(D.5) 

if I is odd and by 



1 

2 

f{6) = -kQl{x) ^{AmOix - Xm) + BmOix + Xm)) , X = COS 9 

(D.6) 

2 

J2{Am + B^) = l , 1 = 2,4,... 

m=l 



if / is even. The A^'s and the -Bm's are arbitrary constants. The constraints on their sums 
to be unity is necessary because at x ~ 1 the solution should behave like / = —kQi{x) ^ 
^ ln(l — x) in order to have the correct normalization. Notice that the solutions (p.5[ ) ( PX6| ) 
are discontinuous at the points Xm and that for Z = there is no solution (unless k = 0) 
because there exist no point where the slope of Qo{x) is zero. It should be noted that we 
can add in ( p.5| ) ( p.6| ) any solution of the homogeneous equation in ( p.l| ) , namely a linear 
combination of the spherical harmonics Yj:^{9, (p), m = —I, — / + 1, . . . , /. 
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Let us next consider the differential equation ( p.l| ) but with Laplacian 

A(2) =a,(x2-l)a, + ^^ , x = coshx, (D.7) 
^ ' — 1 

on the hyperboloid with the Lobachevsky metric ds'^ = dx^ + sinh^ % d(p'^. We consider 

the case when c = l{l + 1), / = 0, 1, . . .. Then the solution is given by 

fix) = kQiix), x = coshx, / = 0,1,.... (D.8) 

Notice that in this case there is a solution for / = 0. Also we did not need to use any 
9-functions because /(x) as given by ( p.8| ) becomes asymptotically zero for large x it 
should be. 

Finally let us present the solution of the equation for the shift function (|3.8| ) which is 
valid in any number of dimensions d > 4. This can be written down using the generalized 
spherical functions ||2^, which are representations functions of SO{d — 1). If the number 



of spacetime dimensions is even they are defined in terms of usual Legendre polynomials 

as 

p/^'^^(cosg) = ^1^^^'^ {^^y-' Pi+,_,{cose) , d^2r + 2, r = l,2... , (D.9) 
whereas if it is odd in terms of Chebishev polynomials as 

{2Tvy ^ a cos 6* smu 

(D.IO) 

The only properties of the generalized spherical harmonics we need here are 

A(,_2)P/'-') = -/(/ + rf-3)P/'-') 

1=0 

Using the above formulae we find that the solution for the shift function is 

Obviously the above solution is not valid if a coincides with one of the eigenvalues of the 
Laplacian, i.e. if a = —n{n + d — 3), n = 0, 1 . . .. Then one has to use the analogs of the 
Qi^s (see |]2^) but we will not elaborate on that any more here. 
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